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The study of the noise induced effects on the dynamics of a chain molecule crossing a 
potential barrier, in the presence of a metastable state, is presented. A two-dimensional 
stochastic version of the Rouse model for a flexible polymer has been adopted to mimic the 
molecular dynamics and to take into account the interactions between adjacent monomers. 
We obtain a nonmonotonic behavior of the mean first passage time and its standard 
deviation, of the polymer centre of inertia, with the noise intensity. These findings reveal 
a noise induced effect on the mean crossing time. The role of the polymer length is also 
investigated. 
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I. INTRODUCTION 



Many important biological processes are gov- 
erned by the transport of molecules across mem- 
branes. In cellular systems proteins translocate 
from the cytosol into the endoplasmatic retic- 
ulum [Hu et al., 2005; Chiam et al, 2006], or 
into mithocondria or chloroplasts [Frappat et al, 
2003]; RNA molecules are transported across a 
nuclear membrane after their synthesis [Schatz 
& Dobberstein, 1996]; in infections from bacteria 
resistent to antibiotics, alternative therapy uses 
bacteriophages which push a long chain of DNA 
molecules through the pores of a membrane [Ad- 
hya &; Merril, 2006]. The study of the dynam- 
ics of translocation of individual polymers across 
nanometer-scale pores is essential to help the un- 
derstanding of how biological systems work and 
it also represents a powerful instrument to in- 
vestigate the structural characteristics and con- 
formational changes of chain molecules, their in- 
teractions with the walls of the pore and spe- 
cific folding or unfolding processes related to the 
crossing event. A fundamental work on poly- 
mer translocation was carried out by Kasianow- 
icz et al. [1996]. In their experiments, single- 
stranded DNA molecules are forced by an ap- 
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plied electric field to move through a 2.6-nm di- 
ameter ion channel in a voltage-biased lipid bi- 
layer; during the passage across the pore, the 
polymer significantly reduces the single-channel 
ion current. By measuring the duration time 
of each current drop - blockade -, they have 
demonstrated that the crossing time linearly in- 
creases with the polymer length. By using a sim- 
ilar experimental setup, Meller et al. [2002] have 
also shown that the time-length relationship de- 
pends on the number of adenines in the polymer, 
concluding that the translocation time is heav- 
ily influenced by specific interactions between 
the DNA and the nanopore walls. They also 
suggest stronger interactions at lower tempera- 
tures because of the longer measured transloca- 
tion times [Meller et al., 2000]. On the other 
hand, very interesting experimental results by 
Han et al. [1999] show longer crossing times for 
shorter DNA molecules, suggesting the existence 
of a quasi-equilibrium state of the polymer after 
being trapped for a time. Theoretical studies 
have demonstrated that the translocation time 
of a long N-segments chain molecule, diffusing 
across a membrane, scales as N^ if there is no 
free energy bias between the two sides, or as N^ 
in case of adsorption [Sung & Park, 1996; Park 
& Sung, 1998]. Strong interactions between the 
polymer and the pore modify the dynamics of 
the portion of the chain inside the channel, mak- 



ing linear the relationship between the crossing 
time and N [Lubensky & Nelson, 1999]. In the 
presence of an externally applied potential, the 
linear dependence of the crossing time on the 
length of the polymer is confirmed only for the 
case of an hairpin crossing mechanism of translo- 
cation [Schatz & Dobberstein, 1996; Sebastian 
& Paul, 2000]. The dependence of the crossing 
time on the length of the polymer chain, as re- 
ported in different experimental works, may be 
influenced by the different geometrical and phys- 
ical characteristics of the pore-channel device, 
the adopted polynucleotides, the intensity of the 
driving electric fields. We are convinced that, in 
same cases, a special condition for the polymer 
travelling across a nanometer-scale channel can 
be realized to make it temporarily trapped into 
a metastable state, as already suggested by Han 
et al. [1999]. For this reason, we decided to start 
a theoretical study of the dynamics of polymers 
in a state out of equilibrium under the influence 
of a noisy environment. In this paper, we firstly 
describe the model and the mathematical frame- 
work adopted to simulate the polymer dynamics 
and, finally, we comment our results. 

II. THE MODEL 

Stochastic modelling of polymer dynamics 
starts with the choice of a model for the chain 
molecule. In the simplest model of a polymer, 
monomers are freely-jointed rigid rods of a fixed 
length (Kuhn length) and any kind of interac- 
tions among them is neglected. For our studies, 
we have adopted the Rouse model for a flexible 
polymer [Rouse, 1953; Doi & Edwards, 1986]. 
This model describes the molecule as a chain 
of N monomers connected by harmonic springs 
(Fig. EI). 




FIG. 1 Schematic representation of the polymer 
chain. 

At this stage, we have neglected the excluded- 



volume effect and any hydrodynamic interac- 
tion between monomers and solvent, as mod- 
elled by the Zimm theory [Zimm, 1956]. In a re- 
cent experiment, a fluorescent correlation spec- 
troscopy technique has been applied to monitor 
the stochastic motion of individual monomers 
within isolated DNA molecules [Shusterman et 
al., 2004]. By comparing the observed time de- 
pendence of the mean-squared diplacement of 
the fluorescent monomer with that predicted by 
the available models, Shusterman and coworkers 
assert that double stranded DNA (dsDNA) be- 
haves like a semiflexible coil following the Rouse 
model, while in the case of single stranded DNA 
(ssDNA) the monomers show a Zimm-type ki- 
netics. A similar analysis of the short-time dy- 
namics demonstrates that intramolecular hydro- 
dynamic interactions in dsDNA chains cannot 
be ruled out [Tothova et al., 2005; Petrov et 
al., 2006], suggesting the application of a mixed 
Rouse-Zimm model [Lisy et al., 2005]. We have 
decided to start our investigation by using a rela- 
tively simple bead-spring model in order to avoid 
any secondary effect of metastability induced by 
the combination of an excessive number of pa- 
rameters. 

The translocation of a long molecule through a 
narrow pore is caused by a forced process of un- 
folding by which the coil is stretched across the 
channel, returning in its native folded state again 
at the exit only. Within the pore region, the free 
energy per segment of the polymer is higher, rep- 
resenting a barrier for the motion of the molecule 
[Sung & Park, 1996]. In the present work, we 
study the effect of the temperature fluctuations 
on the dynamics of a chain polymer escaping 
from a metastable state through a potential bar- 
rier in a two-dimensional domain. The polymer 
motion is modeled as a stochastic process of dif- 
fusion in the presence of a time-independent po- 
tential having a profile described by the following 
equation: 

U{x) = ax^ - 6x^ (1) 

with parameters a=0.3 and 6=0.2, as already 
adopted by Fiasconaro et al. [2003] for the unidi- 
mensional single-particle case. The correspond- 
ing 2D-surface is shown in Fig. [2l The drift of 
the i-th monomer of the chain is described by the 
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following coupled Langevin equations [Gardiner, 
1993]: 



X 



dx _ dU{ 
dt dx 



+ Fint{x) + VD^^{t), 



(2) 
(3) 



where Fint{x) and Fint{y) are the x and y 
components, respectively, of the harmonic 
forces between adjacent monomers and ^x(i) 
and Cy(t) a-re white Gaussian noise modeling 
the temperature fluctuations, with the usual 
statistical properties, namely {^i(t)) = and 
{Ut)^j{t + r)) = 6,j{T) for = x,y). We 
have numerically solved the above system of 
equations by performing a set of 10^ simulations 
for 22 different values of the noise intensity D 
and 5 different polymer length configurations. 
In this study, we put the starting distance Lq 
between two adjacent monomers of the molecule 
chain equal to the rest length of the ideal spring 
connecting them. We have chosen Lq = 0.01 and 
the harmonic constant equal to 1.0 (in arbitrary 
unit). The number of monomers changes 
between 2 and 10. The initial position is set to 
a;o=1.3 for all monomers; every simulation stops 
when the x coordinate of the center of mass of 
the chain reaches the final position at Xf=2.2. 
Then, the Mean First Passage Time (MFPT) 
and its standard deviation are calculated. 





FIG. 2 Potential surface of the two-dimensional bar- 
rier. Inset: projection of the potential surface on the 
plane y = 0. 



III. RESULTS 

First results of this study suggest the pres- 
ence of the noise enhanced stability (NES) effect 
[Mantegna & Spagnolo, 1996; Agudov & Spag- 
nolo, 2001; Fiasconaro et al., 2003; Fiasconaro 
et al, 2005] for a short polymer escaping from 
a potential barrier. In the left panel of Fig. [3] 
we show how the MFPT changes with the noise 
intensity at different polymer lengths. Assum- 
ing identical unstable initial conditions we find 
that: (i) the NES effect is more relevant for 
shorter chains (lower number of monomers); (ii) 
the maximum of the MFPT rapidly decreases at 
increasing polymer length and its position shifts 
towards higher levels of noise. 
The dynamics of shorter chains is much more 
similar to that of a single particle escaping from 
the potential well. As a consequence shorter 
chains show more pronounced NES effect, with 
higher maximum values of the MFPT. Moreover, 
the shift of the position of this maximum is di- 
rectly connected with the polymer length. In 
fact, longer polymers need higher noise inten- 
sity values to move the centre of mass back into 
the metastable state, starting from initial un- 
stable states. Longer chain polymers experience 
more than shorter ones the effect of the poten- 
tial slope in the unstable region. It is important 
to note that the behavior of MFPT at very low 
noise intensities {D < 5 -10"^, in Fig. [3^) should 
have divergent trend because the initial unstable 
states used in our simulations are in the diver- 
gent dynamical regime (detailed analysis on the 
divergent dynamical regime can be found in Fi- 
asconaro et al. [2005]). To find this behavior we 
need very long computer simulation times, which 
become greater as the polymer length increases. 
In fact the probability that the entire sequence 
of A'^ monomers experience the potential well is 
much lower than the single particle, and so more 
difficult to see in a finite time computer simula- 
tion. 

The plot of the standard deviation of the first 
passage time as a function of the noise inten- 
sity D is shown in Fig. [3)3. An increasing of the 
standard deviation means a spread of the pas- 
sage time distribution and this represents a di- 
rect evidence of the noise-induced capture effect 
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FIG. 3 (a): Mean first passage time vs. noise intensity for five different numbers N of monomers in the 
polymer. The starting distance between adjacent monomers coincides with the rest distance Lq — 0.01. 
The initial condition xq = 1.3 is the same for all the monomers. The dashed line indicate the shift of the 
maximmn of the MFPT. (b): Standard deviation of the passage time vs. noise intensity for different polymer 
configurations. 
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of the barrier. In Fig. |3]we report the probabil- 
ity density function (PDF) of the passage time of 
the polymer center of mass for the case of N=4 
and four values of the noise intensity, namely 
D = 0.01, 0.1, 1.0 and 10.0. The PDFs obtained 
with lower noise intensities are characterized by 
a strong asymmetry and high long tails. In this 
figure we can clearly see that, while the most 
probable value follows a monotonic decreasing 
trend with increasing noise intensity, the MFPT 
is strongly affected by the rare events of trapped 
particles giving the long tail of the distribution 
at high values of the passage time (see inset in 
Fig. H]). The shape of the tail significantly influ- 
ences the mean crossing time by increasing its 
value with respect to the mode. For the case of 
N=10, the plot of MFPT vs. D does not show 
any evidence of a NFS effect, but not the stan- 
dard deviation, which shows a non-monotonic 
behavior. 

The plot of the MFPT as a function of the poly- 
mer length (Fig. [5^) at different noise intensities 
shows that: (i) the quasi-deterministic MFPT 
(noise level D = 0.01) does not depend on N; 
(ii) for D > 0.01, the dependence of MFPT on 
the polymer length shows a non-monotonic trend 
with the noise intensity: differences of MFPT 
are very sensitive to the number of monomers 
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FIG. 4 Probabihty density function (PDF) of the 
passage time of the polymer center of mass with N=4 
at four different values of the noise intensity. The 
inserted plot shows the tail at high passage times of 
the PDF for D 0.1. 



when the noise level is low, while they are re- 
duced at higher noise intensities. In Fig. O) it 
is interesting to note the capture effect also for 
the polymer with N=4 and the low level of noise 
D = 0.03. 




2 4 6 810 12 2 4 6 81012 

N N 

FIG. 5 (a): Dependence of the MFPT on the number N of monomers at different noise intensities. Initial 
conditions are the same as in Fig. [31 Red hnes are used to discriminate the ascending trend of MFPTs with 
the decreasing of the noise intensity with respect to the opposite behavior at extremely low noise intesity 
(blue lines), (b): Standard deviation of the passage time vs. N for different levels of noise. 



IV. CONCLUSIONS 

Experimental and theoretical studies on the 
translocation time, taken by a chain molecule to 
cross a membrane, show a strong dependence of 
this time with the length of the polymer. In this 
work we propose a possible reason of this depen- 
dence by investigating the stochastic dynamics 
of a polymer, with different lengths, surmount- 
ing a potential barrier. After been trapped for a 
while, the molecule leaves the top of the barrier 
exhibiting an enhancement of the mean crossing 
time as a function of the temperature fluctua- 
tions. This is a typical signature of the pres- 
ence of a metastable state in the system inves- 
tigated. Our analysis shows a clear effect of 
noise-enhanced stability on the mean first pas- 
sage time of short polymers for initial unstable 
conditions. Our simulations demonstrate that 
longer polymers can travel the same channel dis- 
tance faster than shorter chains. This result is in 
agreement with the experimental measurements 
of DNA mobility carried out by Han et al. [1999], 
but a deeper analysis on the role of initial con- 
ditions on the translocation time is necessary. 
Our findings are also in accordance with the ex- 
perimental results of Meller et al. [2000], that is 
the longer translocation times measured at lower 
temperatures, as we have found near the maxima 



of MFPTs (see Fig.[3i). 
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